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A perfect number is a positive integer n such that n equals the sum of all positive integer divisors 
of n that are less than n. That is, although n is a divisor of n, n is excluded from this sum. Thus 
b — 1+2 + 3 is perfect, but 12 1+2 + 3+ 4 + 6 is not perfect. An ACL2 theory of perfect 
numbers is developed and used to prove, in ACL2(r), this bit of mathematical folklore; Even if 
there are infinitely many perfect numbers, the series below, of the reciprocals of all perfect numbers, 
converges. 
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1 Perfect Numbers 

The smallest perfect numbers are 6 = 2-3 = 2^(2^ — 1), 28 = 4-7 = 2^(2^ — 1), 496 = 16-31 = 2^(2^ — 1), 
8128 = 64 • 127 = 2^(2^ — 1). In each of these examples, the second factor, 3,7,31,127, of the form 
2* — 1, is a prime. The Greek Euclid proved ID page 3]: 

Theorem 1 If 2^ — I is prime, then n = 2^^' (2* — 1) is perfect. 

Primes of the form 2^ — 1 are called Mersenne primes. Thus every new Mersenne prime leads to a 
new perfect number. According to Wikipedia ||5l, less than 50 Mersenne primes are known. The largest 
known Mersenne prime is making 2^^’^^5460|'257,885,i6i _ largest known perfect 

number, with over 34 million digits. It is not known if there are infinitely many Mersenne primes, nor if 
there are infinitely many perfect numbers. 

All perfect numbers built from Mersenne primes are even. The Swiss Euler proved every even perfect 
number is built from some Mersenne prime l|2j page 10]: 

Theorem 1 If n is an even perfect number, then n = 2^ ^ (2^ ~ 1 )> where 2^ — 1 is prime. 

It is not known if there are any odd perfect numbers, but Euler also proved 161 page 250]: 

Theorem ?> If n is an odd perfect number, then n = where p is prime and i, p, m are odd. 

ACE2 is used to verify each of these three theorems. 

If there are only finitely many perfect numbers, then clearly the series 

^ I 

perfect/! 

converges. ACE2(r) is used to verify that even if there are infinitely many perfect numbers, the series 
converges. 
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2 The ACL2 Theory 

In number theory, for positive integer n, o{n) denotes the sum of all (including n) positive integer divisors 
of n. The function o{n) has many useful properties, so the definition of a perfect number is reformulated 
in terms of a lU pages 8-9]: 

perfect(n) if and only if o{n) = 2n 


These six properties of a are among those formulated and proved in ACL2: 

1. p is prime if and only if o{p) = p + 1. 

2. If p is prime, then a(p*) = Y!i=QP' = NNN- 

3. If p and q are different primes, then a(p • q) = cj(p) ■ cj(^). 

4. o{k-n) < o{k) ■ o{n) 

5. If gcd(/:,?i) = 1, then o{k ■ n) = o{k) ■o{n). 

6. If p is prime, then gcd(p^, a(p^)) = 1. 

If n = 2'(2'+^ — 1) is an even perfect number, then the exponent i is computed by an ACL2 term, 
(cdr (odd-2"'i n)), that returns the largest value of i such that 2' divides n. 

If n = p'm^ is an odd perfect number, then p,i,m are respectively computed by the ACL2 terms 

• (car (find-pair-with-odd-cdr 

(prime-power-factors n))) 


• (cdr (find-pair-with-odd-cdr 

(prime-power-factors n))) 


• (product-pair-lst (pairlis$ (strip-cars 

(remove1-equal 
(find-pair-with-odd-cdr 
(prime-power-factors n)) 
(prime-power-factors n))) 
(map-nbr-product 
1/2 

(strip-cdrs 
(remove1-equal 
(find-pair-with-odd-cdr 
(prime-power-factors n)) 
(prime-power-factors n)))))) 


These terms implement the following computation: 

1. Factor n = n5=o ’^he product of powers of distinct odd primes. 

2. Exactly one of the exponents, say eo, will be odd and all the other exponents will be even. 

3. p is the prime with the odd exponent and i is the unique odd exponent. Son = p' ■ 05=1 pY^ ■ 
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4. Then m = 05=1 pf n = p'rrP'. 

ACL2 is used to verify a result of B. Hornfeck, that different odd perfect numbers, ni = p\m\ / 
n 2 = P 2 ^\ have distinct mi ||6l page 251]: 

Theorem 4 Ifni = p'\^\ cind n 2 = P 2^2 perfect numbers and m\ = m 2 , then ni = ?i 2 - 

Theorems 2, 3, and 4 are enough to prove the folklore that the series, of the reciprocals of all perfect 
numbers, converges. 

3 ACL2(r) 

ACL2(r) ||3l is based on Nonstandard Analysis 111 01 which provides rigorous foundations for reasoning 
about real, complex, infinitesimal, and infinite quantities. There are two versions of the reals 

1. The Standard Reals, ^’■R, is the unique complete ordered field. This means fhaf 

• Every nonempfy subsef of fhaf is bounded above has a least upper bound. 

There are no non-zero infinitesimal elements, nor are there are any infinite elements in ^*'R. 

2. The HyperReals, *R, is a proper field extension of ^*'R: ^*'R ^ *R. There are non-zero infinitesimal 
elements and also infinite elements in *R. 

Here are some technical definitions. 

• X G *R is infinitesimal; For all positive r G ^^R, (|x| < r). 

0 is the only infinitesimal in ^'■R. 

(i-small x) in ACL2(r). 

• X G *R is finite; For some r G ^'■R, (|x| < r). 

(i-limited x) in ACF2(r). 

• X G *R is infinite; For all r G ^’'R, (|x| > r). 

(i-large x) in ACF2(r) 

• x,y G *R are infinitely close, x ss y; x — y is infinitesimal. 

(i-close X y) in ACF2(r). 

• Uoo is an infinite positive integer constant. 

(i-large-integer) in ACF2(r). 

Every (partial) function / ; ^’■R” 1 —^'■R^ has an extension */ ; *R" 1 —*R^ such that 

1. Forxi,--- ,x„ G ^^R, A«) = /(^b-- - An)- 

2. Every first-order statement about / true in ^’■R is true about */ in *R. 

Example. 

(Vx) [sin^ (x) -h cos^ (x) = 1] is true in ^’■R. 

(Vx)[*sin^(x) -|-*cos^(x) = 1] is true in *R. 

Any (partial) function / ; ^'■R” 1 —)• ^'■R^ is said to be classical. 

• Identify a classical / with its extension */. 

That is, use / for both the original classical function / and its extension */. 
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• Use (V^*';c) for (V;r G i.e. “for all standard x.” 

Use for (3x G i.e. “there is some standard a:.” 

• “(V.r)[sin^(A:) +cos^(.r) = 1] is true in ^’■R” becomes “(V^’'.r)[sin^(A:) +cos^(x) = 1] is true in *R.” 

“(V.r)[*sin^(.r) + *cos^(.r) = 1] is true in *R” becomes “(V.r)[sin^(.r) +cos^(.r) = 1] is true in *R.” 

Numeric constants, c, are viewed as 0-ary functions, c : ^'■R° i—)> ^’■R or c : *R° i—*R. Thus, elements 
of ^’■R, such as 2,4, —1, are classical. But elements of *R — ^'•R, such as the infinite positive integer Uoo, 
are not classical. Functions defined using fhe nonsfandard concepfs of infinifesimal, finife, infinife, and 
infinifely close are nof classical. 

Let / be a (partial) unary function, whose domain includes the nonnegative integers, into the reals. 
Here are three possible definitions for the real series YT=ofi^) converges. The first two are versions 
of Weierstrass’ traditional definition that the real series converges. One version for the standard reals, 
another version for the hyperreals. 

1. (defun-sk 

Series-Converges-Traditional-Standard ( ) 

> 0)(3^'' integer M > 0)(V^'' integer n)(n > M | — L| < e) 

) 

2. (defun-sk 

Series-Converges-Traditional-Hyper ( ) 

(3L)(Ve > 0)(3 integer M > 0)(V integer n)(n > M ^ \ < ^) 

) 

3. (defun-sk 

Series-Converges-Infinitesimal ( ) 

(3^''L)(Vinfinite integer n > 0 ){Y 4 =of{i) ~ L) 

) 

For classical /, ACL2(r) verifies fhese fhree definifions are equivalenf. ACL2(r) also verifies for classical 
/, wifh nonnegative range, these definitions are equivalent to this nonstandard definition |]T1 : 

• (defun 

Series-Converges-Nonstandard ( ) 

ir=o/(0 is finite 

) 

Recall that the upper limit, n^o, on this is nn infinite positive integer constant. 


4 The Series Converges 


Use the definition, Series-Converges-Nonstandard, to verify, in ACL2(r), the convergence of 
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by showing this sum is finite: 
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Recall rioo is an infinite positive integer constant. 

Verify the previous sum is finite by showing both of the summands on the right side below are finite. 
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By Theorem 2, even perfect numbers, k, have the form k = 2'(2'+^ — 1). Since 2'(2'+^ — 1) > 2', 
2 '( 2 '+'-i) — h- Iriduction on n verifies Ya=o = 2 — Thus for any positive integer, n, including 

n = nooi 
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By Theorem 3, odd perfect numbers, k, have the form k = p‘m^. Since p'rrP' > m^, < By 

Theorem 4, no square, m^, appears more than once in 



perfect(A^) 
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Induction on n verifies Lm-i ^ <2-2 , Thus for any positive integer, n, including n = Ucoi 
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Therefore, for any positive integer, n, including n = nod 
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and 

- is finite. 
k 
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The heart of this proof is that the partial sums 
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are bounded above (by 4). This can be stated and carried out entirely in ACL2. The Reals and ACL2(r) 
are required to formally state and prove the series converges. 
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A ACL2(r) Books 

A.l prime-fac.lisp 

Unique Prime Factorization Theorem for Positive Integers. 

An ACL2 book as well as an ACL2(r) book. 

A.2 perfect.lisp 

Perfect Positive Integers. 

An ACL2 book as well as an ACL2(r) book. 

Over 500 events, incrementally built Summer 2013 - Spring 2015. 

A.3 seriesl.lisp 

The CLASSICAL series, Seri, converges (to a STANDARD real L). 

A.4 seriesla.lisp 

The CLASSICAL NONNEGATIVE series, Seri a, converges (to a STANDARD real E). 
A.5 sumlist-1.lisp 

Some nice events from sumlist.lisp plus additional events. 

A.6 sum-recip-e-perfect.lisp 

The sum of the RECIPROCAES of the EVEN PEREECT positive integers converges. 
A.7 sum-recip-o-perfect.lisp 

The sum of the RECIPROCAES of the ODD PEREECT positive integers converges. 
A.8 sum-recip-perfect.lisp 

The sum of the RECIPROCAES of the PEREECT positive integers converges. 
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